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ABSTRACT 

A massive (^10 2 -10 4 M Q ) black hole at the center of a globular cluster will tidally disrupt binaries that pass 
sufficiently close. Such an encounter results in the capture of one binary component by the black hole and 
ejection of the other star. A 10 3 M Q black hole disrupts one binary every l-10Myr; captured stars orbit with 
initial periods of days to decades. Although >100 stars are captured in only 1 Gyr, the number with periods of 
< 10 yr at any one time is probably less than 10, regulated by stellar interactions and tidal disruption by the black 
hole. Binary disruption and stellar capture should be critically examined with detailed numerical simulations, 
in order to corroborate the analytic results given here and to weigh possible observational consequences. 
Subject headings: black hole physics — globular clusters: general — stellar dynamics 



1. INTRODUCTION 

Black holes of mass ^10 2 -10 4 M Q occupy a special niche 
in several subfields of theoretical astrophysics. Such ob- 
jects might naturall y form as remnants of the first stars in 
th e Universe jB arkana & Loebll200!l iMadau & Reel 1200 It 

f:hneider et al. 2002) , or be built th rough merging of stars 
'ortedes Zwart & McMillan 2002) or stellar-mass black 
holes (Miller & Hamilton 2002) in young star clusters. These 
black holes may play critical roles in the formation and 
gro wth of su permass ive (>10 6 M^) b lack holes in galactic nu- 
clei (Ebisuzaki et a ll2001lllslam. Tavlor. & SilU2004l) . 

Intermediate-mass black holes (IMBHs) have been in- 
voked to resolve several recent observational quandaries. 
Some of the "ultraluminous" X-ray sources detected in other 
galaxies may be accreting IMBHs (e.g, Miller & Colbert 
2004), although t he ma j ority are likely stellar-mass black 
holes l lKing et alJ 1200 ll iRaupaPOrt. Podsia dlowski. & Pfahll 
2004). There have been recent claims of possi- 
ble IMBH detections in two globular clusters (Ml 5 in 
the Milky Way and Gl in M31) based on HST data 
( Gebhardt, Rich . & Hol2002t iGerssen et aT]l2002l). However 
the data do not require IMBHs, as shown by Baumgardt et al. 
( 2003 albl) and iMcNamara. Harrison. & Andersonl J2003I) . 
Colpi, Possenti, & Gualandris (2002) cite evidence for a sin- 
gle or binary IMBH in the cluster NGC 6752 based on the 
anomalous positions and accelerations of cluster radio pul- 
sars. Finally, a single or binary IMBH could account for the 
net rotation observed in the center of M15 (iGebhardt et alJ 
2000; Gersse nendl2002lMiiler & Colberl2004) . 

Because of the broad scientific scope sampled by IMBHs, 
it is crucial to investigate what unique dynamical signatures 
they imprint on their environments. The high stellar densi- 
ties and proximities of Galactic globular clusters make them 
prime locales to search for IMBHs, as emphasized by much of 
the recent observational work; for a review of ear ly theoret- 
ical and observational studies, see Shapiro ( 1985). Here we 
introduce a process that deposits stars directly into compact 
orbits (e.g., <100 AU) around a massive cluster black hole. 

Bina ry stars are abundant in globular clusters (e.g., 
iHut et alJl99l . A binary that passes sufficiently near a cen- 
tral black hole is disrupted by tidal forces, resulting in the cap- 
ture of one star and ejection of the other (see §|2jl. Thi s type 
of inte racti on was studied in the context of Sgr A* bv iHillsl 
( 1988) and lYu & Tremainel ( 120031) as a mechanism for accel- 



erating stars to very high speeds, and by iGould &"O uillen 
(120031) to account for the young stars near Sgr A* (e.g., 
Eisenhaue r et alJl2005l) . The disruption rate in globular clus- 
ters is estimated in § [3] and the fate of captured stars is dis- 
cussed in § @] Future directions for studying this mechanism 
in more detail are suggested in §|5j as are potentially observ- 
able consequences. 

2. BINARY DISRUPTION 

Consider a binary system of mass Mb = M\ +M 2 with semi- 
major axis cih on a parabolic approach to a black hole of mass 
M. ^> Mb- Let r be the distance between the binary barycenter 
and the black hole, and s be the minimum distance. The dif- 
ferential acceleration across the binary (~GM,fli,/r 3 ) exceeds 
the binary self-gravity (~GMj/flj) within the tidal radius, 



r,=a b &f m /M b ) 1/3 = 10a b Qj 



>V3 



(1) 



where M m /Mb = 10 3 (23- For a disrupted binary, the tidal per- 
turbation is essentially impulsive, and the additional relative 
velocity imparted to the binary components has a magnitude 
of ~Vj, = (GMb/ab) 1 ' 2 . Immediately following the disruption, 
the two stars are mutually unbound and move relative to the 
black hole with the original systemic velocity of magnitude 
u, ~ (2GM./r,) 1//2 , plus contributions from the prior orbital 
motion and tidal perturbation. Random orbital orientations 
and phases at the time of disruption lead to a range of orbital 
parameters and asymptotic speeds for stars that are captured 
or ejected; below are estimates for the typical values. 

If star 1 is captured, we can estimate the expectation value 
of its specific energy, E\, relative to the black hole by as- 
suming that its resultant velocity after the tidal perturbation 
is directed opposite to the binary systemic velocity and has a 
magnitude of (Mi/Mb)v b - To leading order, we obtain 
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The corresponding semimajor axis and orbital period are 



a x =-GM./2E x ~ 35a h (M h /M 2 )Q 



2/3 
3 ! 
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and 



P l ~lP b {M b /M 1 ? l2 Q 



1/2 



(4) 



where P b is the binary period. The pericenter separation of 
star 1, a\(\-e\), must be comparable to s, so that 



which is >0.7 for £>3 > 1. Star 2 has a typical specific energy 
of E% — (Mi/Mb)UtVb, and thus its speed at infinity is 



v 2 ~Av h (M l /M b yi 1 Q 



1/2^1/6 



(6) 



which may be >100kms 1 for typical binary parameters. A 
star ejected by this process can easily escape its host cluster. 

3. TOTAL DISRUPTION RATE 

We now employ "loss-cone" theory to estimate the rate of 
binary disruptions. The approach taken here is not very rigor- 
ous, but captures the essential physics. More pr ecise accounts 
of loss-cone dynamics can be fou n d in, e.g.. iFrank & Reesl 
lH^7^.lLightman & ShaDirollfT977l).ICohn & Kulsradl ill 9781) . 
ISver & IJlmeil ill 999ft . and Magorrian & TremaineT Il 999ft . 

A binary with semimajor axis a b at radius r 3> r, is tidally 
disrupted if its specific angular momentum about the black 
hole is less than J, = {2GM t r t ) x l 2 . If the local one-dimensional 
stellar velocity dispersion, er(r), is adopted as a typical binary 
systemic speed, the disruption condition defines a double- 
sided loss cone in velocity space with semiangle 6 t = J t jra. 
Perturbations by distant stars cause the binary barycenter to 
undergo a random walk in phase space. In a dynamical time, 
Td i~ r/a, the velocity is deflected by an angle 6^ ~ (tj / r, ) 1 / 2 , 
measured at radi us r. The local two-b ody relaxation timescale 
is given by (e.g., Spitzer & H aril971l) 



Tr(r) = 



o\r) 



(cr/lOkms" 1 ) 3 
50 . ' - — - Myi. 



AG 2 m 2 n(r) Ajo(n/10 5 pc- 3 ) ' 



(7) 



where n(r) is the stellar number density, A = A10 10 is pro- 
portional to the Coulomb logarithm, and we have assumed 
m = 1 Mq for the mean stellar and binary mass, which is fixed 
in the subsequent calculations. 

The ratio 8^ / 6, increases monotonically with r for profiles 
n(r) and a(r) suitable to globular clusters. Two loss-cone 
regimes are loosely demarcated by 9d/0\ c = 1, defining a crit- 
ical radius r c . For r > r c , binaries can scatter into and out of 
the loss cone in a dynamical time, and the loss cone is full sta- 
tistically. If r < r c , binaries within the loss cone are disrupted 
in a dynamical time, and the disruption rate is limited by the 
diffusion of binaries across the loss-cone boundary. 

We assume that the cluster center consists of a power-law 
density cusp around the black hole, embedded in an isother- 
mal cluster core of radius ro. The cusp is inside the radius of 
influence of the central black hole: 



n = GM./a 2 ~ 0.04M 3 er^pc 



(8) 



where M, = IQ^M^Mq, and 07 = lOeriokms 1 is the velocity 
dispersion at r,. Each zone has different power-law profiles, 



n(r) = rii(r/ri) a and a 2 (r) = a Hr/r,) ° , where a = 8 = in 
the core, and a ~ 1.5-2 (e.g., Ba hcall & Wolll977t lYoung l 
H97l) and/3= 1 in the cusp. The critical radius is 



r c /r i = {a b /a c flW-°'\ 
where, given «,■ = 10 5 «5pc~ 3 , we have 

a c ~ 0.1 n 5 M^ /3 cT^ AiqAU 



(9) 



(10) 



Cluster binaries interact frequently with single stars and 
other binaries. These encounters shape the binary population 
as the cluster and its stars evolve, a complex problem that can 
only be ful ly unrave l ed wi th detailed cluster simulations (e.g., 
Fregeau et al. 2003; Giersz & Spurzem 2003). Here we as- 
sume for simplicity a differential distribution p{a b ) cx a~ b over 
the range a b = 0.01-1 AU, and a fixed binary number fraction 
f b . As the cluster evolves, the shape and limits of p(a b ) vary, 
and f b in the cor e may decrease from >0.5 initially to current 
values of <0.1 dlvanova et al. 120051) . 

When the loss cone is full, a fraction ~df of the binaries 
with semimajor axis a b at radius r are disrupted in a dynamical 
time. The differential rate is 



dUm (M. 
— ~ 4nf b GM. — 

aa b \ m 



1/3 



p(a b )a b 



dr n(r) 
r a{f) 



(11) 



The contribution from r > r, dominates, and we obtain 



Tim - 10" 5 f b n s (a b (AXJ))A4 /3 (Tio In ( - ) yr" 1 , (12) 



where (a b {AU)) is the mean semimajor axis in AU. 

In the diffusive regime, the fraction of binaries disrup ted in 
a time T>(r) is ~l/ln(2/0,) (e.g. JFrank & Reeslll976f) . The 
differential diffusive disruption rate is then 



:4-Kf b G m p(a b ) / , (13) 



da b 



cr 3 (r) 



where the weakly varying factor A/ln(2/0 f ) ~1 has been 
dropped. The lower limit in the above radial integral is rel- 
atively unimportant if it is <Cr c . Numerically, we find 



ftdiff-io-V/^^roV -1 



(14) 



1 for both the 
= 1. The total 



A binary fraction of f b ^0.1 and mean semimajor axis of 
(a b ) ~ 0.1 AU give disruption rates of ~10~ 7 yr" 
full and diffusive regimes when M3 = n$ = crjo 
disruption rate is then ^10~ 7 -10~ 6 yr _1 for the same param- 
eter values. For a b = 0.01-1 AU and M3 ~ 1, captured stars 
have a\ ~ 1-100 AU when M\ ~ M 2 , and periods of ^10d to 
^30yr. From eq. 0, p(a\) should have a similar shape as 
p(a b ). Our choice of p(a b ) cx a~ b x implies that the distributions 
of log flj , log Pi , and log V2 should be roughly flat. 

As an illustration, we computed 10 4 scattering encounters 
using a binary-singl e scatter i ng cod e similar in logic to that 
described in Hut & Bahcall ( 1983). The motion was inte- 
grated with the symplectic algorithm of Mik kola & Tanikawal 
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FIG. 1 . — Simulated outcomes for binaries with Mi = M2 = 1 Mq disrupted 
by a IOOOMq black hole. The panels to the top and right of the scatter plots 
show the histograms of logaj, logV2, and 1— e\. Below the dashed curve 
(a\(\ — ej) = IOKq) a solar-type star is tidally disrupted by the black hole. 



( 119991) and Pre t o & Tremaind(ll999l) . We chose a black-hole 
mass of 1OOOM and stellar masses of M\ = M 2 = 1 M . The 
binary orbits were initially circular, with separations drawn 
from p(at,) oc (ab = 0.01-1 AU). The speed at infinity was 
fixed at Voo = lOkms" 1 . Impact parameters, b, were drawn 
from a distribution that is flat in b 2 ~ 2sG(M m +M b )/v 2 00 , 
limited so that s < 2r,. Figure 1 shows that the distribu- 
tions of logfli and log V2 are roughly flat in the middle, with 
high and low tails not predicted in § [2] Solar-type stars with 
ai(l-ei)< IQRq (below the dashed curve in Fig. 1) are prob- 
ably tidally disrupted by the black hole (see eq. fl|). 

4. DYNAMICAL EVOLUTION OF CAPTURED STARS 

Each binary disruption by a ~10 3 M© black hole deposits 
one star into a compact, eccentric Keplerian orbit. In 100 Myr, 
~10— 100 stars may be captured. However, a collection of 
>10 stars orbiting the black hole with periods of <10yr will 
not persist in the time between captures. Mechanisms that 
drive the evolution of captured stars include gravitational and 
collisional two-body encounters, and tidal disruption by the 



black hole. These processes are quantified below, under the ad 
hoc simplifying assumptions that the stars have equal masses 
of m = IMq, follow a power-law density profile with a ~ 2, 
and everywhere have a nearly isotropic velocity field. We 
neglect relatively weak and rare perturbative encounters be- 
twee n captured stars and cluster stars on highly elongated or- 
bits ( GultekinlMniiEI&HaMto32003). 

Consider N point masses orbiting the black hole. Let r/, 
be the radius that contains half of the stars and Pi, be the 
period of a circular orbit at r;,. The relaxation time at r/, is 
r h - lO^M^Ar 1 (eq. with A = ln(M./m ) ~ 10). "Res- 
onant relaxation" ( Rauch & Tremaine 1996) of a star's or- 
bital angular momentum occurs on a shorter timescale of 
T T h es ~ 10 3 f/,M3, driving the collection of captured stars to- 
ward velocity isotropy without affecting the energy distribu- 
tion. If Ph ~ lOyr and Af > 10, the orbits spread significantly 
in the radial direction in the time between captures. Cap- 
tured stars on the widest orbits will mingle with other (non- 
captured) stars within r,-. The outward transport of captured 
stars into the cusp may noticeably alter the cusp structure. 

A star may be unbound f rom the black hole follo wing a sin- 
gle strong encounter (e.g., Lin & Tremaine 1980). An upper 
limit to the ejection rate is obtained by setting the maximum 
two-body impact parameter to £> e j = Gm/a 2 (r), characteris- 
tic of a ~90° deflection but not necessarily an ejection. The 
local rate per star is ^Tm{r)b 2 ^a(r). Integrating over all cap- 
tured stars, we find a rate of <10~ 7 P^ 1 Mf N 2 . For AT ~ 10 
and P/ t ~ lOyr, most stars may be ejected in the time between 
captures. Under the point-mass assumption, energy conser- 
vation demands that ejection and relaxation leave at least one 
star tightly bound to the black hole. When finite stellar radii 
are considered, two additional processes must be addressed. 

First, the (non-resonant) rate at which the black hole tidally 
disrupts captured stars is ~N/Tf, ~ 10~ 5 P^ 1 M^ 2 N 2 , which is 
much more efficient than ejection at removing stars. Resonant 
relaxation may increase the rate substantially, since angular 
momentum diffusion drives stars into the loss cone. Note 
that if a considerable fraction (e.g., >10%) of the captured 
stars are ultimately disrupted by the hole, the implied rate 
may exceed the swallowing rate for single stars determined 
by previous au thors (~10 yr" 1 for Mi, = n% = crm = 1; e.g., 
Frank & Reeslll976ULightman & Shapiro! 19771) . Second, the 
collision rate for stars of radius R is ^2ttP^ 1 (R/ 'rh) 2 N 2 ~ 
(10- 9 -10- 6 )(^/ J R Q ) 2 A? 2 yr- 1 for r h = 10-50 AU. A collision 
between two main-sequence stars will result in coalescence 
unless the relative speed at infinity is significantly larger than 
the es cape speed from the stellar surface (e.g., Freitag & Benzl 
2004), which may occur within ~1AU of a 1O 3 M black 
hole. 

When N > 10 and P/, ~ lOyr, all of the above processes 
act on timescales shorter than or comparable to the time be- 
tween stellar captures. Thus, at any time, only a few stars are 
expected to orbit the black hole with periods of <10yr. 

5. CONCLUSIONS AND OUTLOOK 

A 10 3 M Q black hole at the center of a globular cluster 
may tidally disrupt thousands of binaries over lOGyr, each 
encounter leaving one star within ^100 AU of the hole. At 
any given time, the number of stars in such compact orbits is 
probably only a few, limited by interactions between captured 
stars and tidal disruption. Before we can assess possible ob- 
servational consequences of these processes, numerical simu- 
lations are required to properly calculate the binary disruption 
rate and the dynamical evolution of captured stars. 
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The clust er evolution codes of [Fregeau et al.l d2003l Monte 
Carlo) and iGiersz & Spurzeml d2003L gas-dvnamical/Monte 
Carlo) are well suited to the problem of determining the bi- 
nary disruption rate. These codes efficiently treat a realistic 
number (>10 5 ) of stars and binaries, and include explicit in- 
tegration of close binary-single and binary-binary encounters. 
The bl ack hole would be rep resented by a central point mass, 
as in lFreitag & Benzl J2002I) . but the imposed spherical sym- 
metry demands that it be fixed at the origin, neglecting its 
Brownian motion. These algorithms cannot adequately fol- 
low the dynamics of a small-A^ system of captured stars. 

Direct A^-body integration has been used to study the 
dynamics of single stars near a m assive black hole 
jBaumgardt et alJ2004bl:lPreto et all2004l) . Such simulations 
are prohibitively time-consuming when N > 10 and there is 
a sizable binary fraction. Simulations with N ~ 10 4 may be of 
value, since the dynamics of captured stars are computed cor- 
rectly. It may also be useful to follow the collisional A^-body 
evolution ~10 captured stars in isolation from the cluster. 

The high binary disruption rate and compact orbits of cap- 
tured stars increase the likelihood that a cluster IMBH will be 
detected, perhaps in one of the following ways: 

1. If, following a stellar tidal disruption, ~0.5M Q is ac- 
creted by the black hole with a 10% mass-to-light conver- 
sion efficiency, the duration of the Eddington-limited flare is 
~10 Af^yr. Such a flare might be visible as a very bright 
UV/X-ray source in one of the ^100 globular clusters in the 
Milky Way if most contain a ~10 3 M Q black hole that dis- 
rupts stars at a rate of ~10~ 6 yr _1 ; but we do not, and both 



suppositions may be highly optimistic. 

2. If the orbit of a captured star with a pericenter 

distance of ~IO.M3.R0 circularizes without the star be- 
ing destroyed, stellar evolution may drive a long (>10 8 yr) 
phase of mass transfer, with accretion r ates of <10~ 8 Afo yr" 1 
and luminosities of <10 38 ergs~ 1 llHopman et al.l 120041 
Portegies Zwart et al. 2004). Steady mass transfer from a 
>10M Q donor has been proposed as a fuel source for ultralu- 
minous X-ray sources in young star clusters. 

3. Hundreds of cluster neutron stars may acquire binary 
companions dynamically, accrete mass, and ultimately turn 
on as millisecond radio pulsars. A ~10 3 M© black hole may 
capture 1 pulsar from a binary system every ~1 Gyr, making 
the detection of such a pulsar unlikely. The timing properties 
of an orbiting pulsar would clearly reveal the black hole. 

4. A captured neutron star or white dwarf that migrates to 
within ^IRq of a 10 3 M Q black hole would produce a de- 
tectable gravitational-wave signal for LISA from distances of 
<100Mpc in th e short time (5, 1 00 yr) before coalesc ence (see 
Baumgard t et all2004btlPortegies Zwart et alJ20"04T) . 
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